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I.  INTRODUCTION 


In  the  Born-Oppenheimer  approximation  for  molecular  collision 
problems  the  electronic  potential  energy  serves  as  the  potential 
energy  for  motion  of  the  molecular  nuclei.  To  facilitate  dynamics 
calculations  for  these  problems  an  analytical  potential  surface  is  a 
very  valuable  tool.  An  analytical  surface  provides  the  necessary 
derivatives  of  the  potential  with  respect  to  the  internuclear  coordi- 
nates over  all  space.  An  analytical  surface  must  also  be  flexible 
enough  to  provide  a reasonable  fit  to  the  known  properties  of  the 
potential  surface  derived  from  experiment  or  ab-initio  calculations. 

Many  molecules  or  molecular  complexes  are  characterized  by  having  their 
minimum  electronic  potential  energy  in  a bent  configuration  (0*0,  180°; 
see  Figure  1).  Such  molecules,  of  which  water  is  an  example,  cannot  be 
modeled  using  a standard  LEPS  (London,  Eyring,  Polanyi,  Sato)  surface. 

In  this  paper  a modified  LEPS  surface  is  presented  which  can  model  the 
bent  configured  water  molecule.  While  not  as  accurate  as  some  other 
1 2 

analytical  methods  ’ the  modified  LEPS  surface  is  easier  to  apply  to 
dynamics  calculations.  Modified  LEPS  surfaces  have  been  developed  for 
other  systems,  HO^  and  CO^,  in  addition  to  water.  Dynamics  calculations 

are  already  complete  for  the  HO^  system. 

II.  MODIFIED  LEPS  SURFACE  FOR  WATER 

An  LEPS  surface  which  approximates  the  equilibrium  ground  state 

for  water  (X^Aj)  has  been  developed.  A single  Sato  parameter  which  is 

a function  of  the  internuclear  coordinates  is  used  to  alter  the  LEPS 
so  that  the  water  molecule  in  and^near  the  equilibrium  configuration 
can  be  modeled.  An  earlier  study  suggested  the  use  of  large  values  of 
the  Sato  parameter  to  produce  wells  or^barriers  along  the  reaction  path. 
The  general  form  of  the  LEPS  potential  is 

v = Qj7  + Q ' + Q ' - [a^2  + a2/2  + a ' 2 


T.  K.S.  Sorbie  and  J. N.  Murrell,  Molecular  Physios  29  (1975),  1287 . 

2.  S.  Farantos,  E.C.  Leisegang,  J.N.  Murrell  and  K.  Sorbie,  Molecular 
Physics  34  (1977)  947. 

3.  A.  Gauss  Jr.,  J.  Chem.  Phys  65  (1976)  4365. 

4.  J.T.  Muckerman,  J.  Chem.  Phys.  56  (1972)  2997. 
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where 


Q>Q./(1*A)  = D?/4(l+A) 

1 1 L 

X { (3  +A)  exp  - 20i(ri-rio) 

- (2  + 6A)  exp  ~Bi(ri-ri0)  } 

and 

^'=0^/(1  + A)  - D®/ 4(1  + A) 

X | (1  + 3 A)  exp  -2&i (ri-ri0) 

- (6  + 2 A)  exp  'Mri.~rio^  ’ 1 _ 1,2> 3’ 

0 

The  Morse  parameters  (D^,  2j>r^g)  are  displayed  in  Table  1.  The 
Sato  parameter  is  given  by 

A = A { f (0) } \ (sin  2 Th)m  ( 

X | (0.94r1r2/r2iQ)Bexp  ^-0.94(r1r2/r21())B  )}, 

where  Th  = Arctan  (rj/r2)  anc^  t*ie  parameters  A,m,  and  B are  adjustable. 

Table  1. 


0-H 


H-H 


r = r 
10  20 


106.53  kcal/mole 
2.29  A"1 
= 0.9706  A 


D®  = 109.5  kcal/mole 
B3  = 1.942  A-1 
r30  0. 7417  A 


This  form  of  the  Sato  parameter  could  be  applied  to  other  three- 
atom  complexes  where  two  of  the  atoms  are  identical.  It  has  the 
proper  symmetry;  that  is,  it  is  invariant  to  interchange  of  r^  and 

r2»  A single  Sato  parameter  has  been  employed  here.  Single  Sato 

parameters  have  been  used  successfully  by  other  workers  on  three-atom 
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4 5 

systems  where  two  of  the  atoms  are  identical  ’ . The  Sato  parameter 
ranges  from  near  zero  for  r^  = r^,  -*■  “(or  r2  = t2q>  rj  “0  to 

slightly  greater  than  two  near  the  minimum  energy  configuration 
(r:  = r2  * 0.87  A,  0 * 106°). 

The  function  f(0)  contains  additional  adjustable  parameters.  In 

6 3 4 

the  case  of  the  water  molecule  f(0)  = [sin(Thl)]  + y [ sin(Th2)J  , where 

6 1 ^ 

Thl  = r 0 and  Th2  = y 0.  The  exponents,  the  constant  multiplying  the 
' 4 t * 

I s in (Th2) ) term,  and  the  constants  multiplying  0 in  Thl  and  Th2  are  the 
adjustable  parameters.  These  parameters  were  varied  so  that  the 
potential  minimum  was  at  0 = 106°,  close  to  the  actual  angle  of  104-105° 
In  addition  f(0)  was  adjusted  so  that  the  difference  between  the  depth 
of  the  potential  well  in  the  linear  configuration  and  the  depth  in  the 
106°  configuration  was  = 1.8  eV  (compare  figs.  2 and  3),  close  to  the 

1.5  eV  difference  shown  in  the  potential  plot  by  Herzberg  ^ . 

The  second  bracketed  term  in  A,  (sin  2Th)  m,  prevents  the  potential 
well  from  extending  to  large  values  of  r^  and  r2  along  the  entrance  or 

exit  valleys;  it  also  helps  eliminate  the  double  minima  near  the  center 
of  the  potential  well  (a  problem  when  large  (>  1)  values  of  the  Sato 
parameter  are  employed).  A large  value  of  the  exponent,  m = 16,  was 
needed  for  the  water  potential  surface. 

The  third  bracketed  term  in  A eliminates  the  large  negative  values 
of  the  potential  that  occur  for  small  values  of  r ^ and  r^  when  A is 

large.  This  term  approaches  zero  as  r^  and  r2  approach  zero  causing 

the  Sato  parameter  to  approach  zero  and  the  potential  (V)  to  form  a 
high  wall  (see  Figures  2,3  and  4).  It  goes  to  zero  for  large  values 
Tj  and  r2  due  to  the  exponential.  By  adjusting  B one  can  ensure  that 

the  potential  minimum  will  be  near  r^  = r2  % r10'  For  The  water 

molecule  B=4  was  chosen.  For  this  1^0  surface  the  minimum  is  near 

r^=r^  * 0.87  A at  0 = 106°;  the  actual  surface  has  its  minimum  at 

rj=r2  0.96  A (at  0 = 104-1050)6 7. 

The  parameter  A multiplying  the  whole  expression  is  chosen  to 
insure  that  the  depth  of  the  potential  well  (at  0 » 106°)  is  close  to 

2 

the  experimental  value.  In  the  case  of  the  water  molecule  A = 10 
was  chosen.  This  value  makes  the  well  depth  * 5 eV  below  the  zero  of 


6.  M.  Baer , J.  Chem.  Phys.  60  (1974)  1067. 

6.  B.J.  Rosenberg , W.C.  Ermler  and  I.  Shavitt , J.  Chem.  Phys.  66 
(1976)  4022. 

7.  G.  Herzberg , Electronic  spectra  of  polyatomic  molecules  (Van 
Nostrand , Princeton 3 1966). 
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(contours  in  eV) 
ep  at  this  angle 


■ 


2 2 

the  potential,  the  separated  OH  ( II)  + H ( S)  state.  Thus  the  overall 
binding  energy  is  * 9.5  eV,  which  is  close  to  the  experimental  estimate 

of  10  eV8. 

In  the  0 + Hy  channel  the  surface  does  not  separate  properly.  It 

l 3 

should  go  to  0(  D)  + H2;  however,  it  actually  separates  to  0(  P)  + H2. 


III.  DISCUSSION 

The  modified  LEPS  surface  for  water  makes  use  of  only  a single 
variable  Sato  parameter.  Since  there  are  two  diatomics  involved,  OH 
and  0.,,  two  Sato  parameters  could  have  been  used  and  would  provide 

additional  flexibility.  Note  that  the  force  constants  for  H->0  were  not 

fit  using  this  surface.  By  using  two  Sato  parameters  matching  the  force 
constants  of  the  LEPS  surface  to  their  true  values  should  be  possible. 

The  two  Sato  parameter  modification  has  already  been  applied  to  HQ2 

with  good  results  for  the  force  constants  of  the  complex. 

The  form  of  the  Sato  parameter  used  here  is  not  the  only  possible 
functional  form.  Other  forms  could  be  adopted  and  will  be  in  the 
future.  In  all  cases,  however,  the  functional  form  must  have  the 
symmetry  of  the  molecule  or  molecular  complex  in  question.  In  summary, 
the  modified  LEPS  has  already  proved  itself  a useful  analytical  potential 
surface  and  with  future  modification  will  be  further  improved. 


8.  B.J.  Rosenberg  and  I.  Shavitt , J.  Chem.  Phys.  63  (1975)  2162. 
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